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ABSTRACT 

We analyze Faraday rotation and depolarization of extragalactic radio point sources in the 
direction of the inner Galactic plane to determine the outer scale and amplitude of the rotation 
measure power spectrum. Structure functions of rotation measure show lower amplitudes than 
expected when extrapolating electron density fluctuations to large scales assuming a Kolmogorov 
spectral index. This implies an outer scale of those fluctuations on the order of a parsec, much 
smaller than commonly assumed. Analysis of partial depolarization of point sources indepen- 
dently indicates a small outer scale of a Kolmogorov power spectrum. In the Galaxy's spiral 
arms, no rotation measure fluctuations on scales above a few parsecs are measured. In the in- 
terarm regions fluctuations on larger scales than in spiral arms are present, and show power 
law behavior with a shallow spectrum. These results suggest that in the spiral arms stellar 
sources such as stellar winds or protostellar outflows dominate the energy injection for the turbu- 
lent energy cascade on parsec scales, while in the interarm regions supernova and super bubble 
explosions are the main sources of energy on scales on the order of 100 parsecs. 



Subject headings: ISM: magnetic fields - 
techniques: polarimetric — turbulence 

Introduction 
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Turbulence in the ionized phase of the inter- 
stellar medium (ISM) of the Milky Way is well 
described on small scales, while its properties on 
larger scales are more uncertain. On scales smaller 
than ~ 10 11 m (~ 10~ 5 pc), the turbulence in the 
ionized medium is well characterized from diffrac- 
tive and dispersive processes in the ISM influenc- 
ing pulsar signals. The comprehensive study by 



Armstrong et all (|1995j ) showed that on scales of 
10 5 to I0 13 m (~ 10" 11 pc to 10~ 3 pc) the power 
spectrum of electron density n e is well described 
by a power law with a spectral index consistent 
with the Kolmogo rov spectral index a — 5/3 
([Kolmogorovl 1 1 94 lh . Most observed electron den- 
sity power spectra are comp atible with a Kol- 
mogo r ov power spectrum (e.g . Spangle r fc Gwirrn 



1990; IStinebring et alj [200 



Wan g et al 



20051 



lYou et al.l 120071 ). although other spectral i ndices 
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have been reporte d (e.g. lLohmer et al.l 12001 



Shis hov et al. 2003). On larger scales the slope 



and extent of the electron density power spectrum 
are much more uncertain. I Armstrong et al. I (|l995l) 
include measurements of H I and of Faraday ro- 
tation measures for which uncertain assumptions 
about the correlation between H I and n e and 
about magnetic fields, respectively, needed to be 
made. Although these results suggest a power 
law that connects to the Kolmogorov power spec- 
trum on small scales, these assumptions make the 
behavior of the power spectrum at larger scales 
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somewhat speculative. 

Fluctuations in the magneto-ionized medium 
on parsec scales have been measured using struc- 
ture function^] of ro tation measures (RMs) o f 
extragalactic sources. Simonetti fc Cordesl ( 1986h 
showed that the structure function of high-latitude 
sources is flat, indicating that fluctuations in 
RM only exist on size scales smaller than the 
scales they probe, viz. about 3°. This means 
that there is no contribution from large-scale fluc- 
tuations from the Milky Way in the RM (al- 
though a constant RM from the Milky Way is 
not ruled out), and that it is most likely the RM 
contribution intrinsic to the sources that domi- 
nates. This was confirmed bv lLeahvl (|l9871 ). RM 
structure functions at lower latitudes, however, 
show structure functions consistent with a power 
law, although the slope te nds to be shallower 
than a Kolmogorov sl ope ( Simonetti fc Cordesl 
1986; Idegg et al.l 119921 ). Other observations in 



or near the Galactic plane ( Sun fe Hani [2004; 
Haverkorn et al. 120031 l2006al ) also find shallow 



slopes for RM structure functions. 

Measurements of the outer scale of fluctua- 
tions, i.e. the scale at which a structure function 
saturates, differ consid e rably across the Galaxy. 
Lazarvan fc Shutenkovl (1990) used the autocor- 



relation function of synchrotron radiation to de- 
rive a typical scale of 90 pc at a distance of 
1 kpc in a region near the north Galactic pole. 
Haverkorn et al.l ( 2006a ) studied sources in the 



Galactic plane and found that the outer scale 
of fluctuations is smaller than about 10 pc for 
the spiral arms, while it is roughly 100 pc for 
the inter arm regions. A similar outer scale of 
~ 90 pc in the magneto-ionized medium is found 
in the Large Magellanic Clo ud using Faraday ro- 
tation (jGaensler et al. 20051) . whereas H I mea- 
surements indicate a much larger outer scale 
of a few kpc in the Large Magellanic Cloud 
(Elmegrc en et al.l 120011). Small Magellanic Cloud 



( Stanimirovic et al 



ies (jWcstpfahl et al 



1999), and in external galax - 



1999; Elmegreen et al.ll2003h . 



7 A structure function measures the amount of fluctuations in 
a quantity as a function of the scale of the fluctuations. The 
second order structure function of a function / is defined as 
D f (56) = {(f(8)-f(6 + 88)) 2 ) g , where 8 is the position ofa 
source in angular coordinates, 89 is the separation between 
sources, i.e. the scale of the measured fluctuation, and ()g 
means the averaging over all positions 8. 



Han et al.l (|2004l ) estimated a magnetic energy 
spectrum with a slope of —0.37 up to scales of 
15 kpc. They suggest that a magnetic energy spec- 
trum which is flatter than Kolmogorov on scales 
larger than the injection scale of 10 — 100 pc is 
dictated by magnetic helicity inversely cascading 
up from the injection scale to larger scales. How- 
ever, the pulsar rotation measure and dispersion 
measure data they use for the power law fit has a 
scatter of several orders of magnitude, making the 
resulting spectrum uncertain. 

Stellar sources of energy input are expected to 
dominate the turbulent driving in the Milky Way, 
except in the outskirts of the Galaxy where star 
formation is low and gravitational sources and 
instabi lities such as the magne t o-rotational insta- 
bility (|Balbus fc Hawleylll99ll: lHawlev fc Balbusl 
1991 ) come into plav (jSelwood fc Balbusl Il999). 
The stellar sources include supernovae, super- 
bubbles, stellar winds, protostellar outflo ws and 
H II regions. iNorman fc Ferrar 2 (|l996l) calcu- 



late a broadband source function mostly dom- 
inated by supernovae, w hich is confirmed by 



iMac Low fc Klessenl(|2004l ). The lNorman fc Ferrar a 
( 19961 ) source function shows a contribution of 
superbubbles to the turbulent driving on scales 
above ~ 100 pc to about a kpc. Turbulent driving 
on these scales is not observed, possibly because 
of the finite extent of the Galaxy So a maximum 
driving scale of about 100 pc due to supernovae 
is often implicitly assumed, although it is not un- 
likely that these sources inject energy into the 
medium (also) on smaller scales. 

What, if any, is the relation of the electron 
density and magnetic field fluctuations at these 
large scales to the Kolmogorov-like spectrum be- 
low 10 13 m? Can both sets of observations be 
reconciled with one spectrum from kilometer to 
parsec scales? What are the characteristics of the 
fluctuations in the ionized ISM on larger scales? 

We address these questions in this paper, using 
RM structure functions from extragalactic sources 
behind the inner Galactic plane. These data are 
discussed in Section [2] We determine the outer 
scale of fluctuations using two independent meth- 
ods: (1) the amplitude and slope of the structure 
function indicate an uncommonly small outer scale 
of Kolmogorov turbulence in the magneto-ionized 
ISM, as we will explain in Section [3j and (2) the 
same conclusions can be drawn from analysis of de- 
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polarization of extragalactic point sources by the 
Galactic ISM, as shown in Section [4j Section [5] 
gives arguments for a steeper (Kolmogorov) power 
spectrum on small scales, and a discussion of the 
results can be found in Section [SI Section [7] pro- 
vides a summary and conclusions. 

2. Data analysis of polarized extragalactic 
point sources 

The data used are from the Southern Galac- 



tic P l ane Survey (SGPS, iMcClure-Griffiths et al 



2005 : Haverkorn et al. 2006f ). a neutral hydrogen 
and full-polarization 1.4 GHz continuum survey of 
the Galactic plane. The continuum part spans an 
area of 253° < I < 357° and \b\ < 1.5° and con- 
tains 148 polarized sources of which the rotation 



meas ure is measured unambiguously (jBrown et al 
20071 ) . The data were obtained with the Australia 



Telescope Compact Array (ATCA) and are pub- 
licly available!!. For more details on the data re- 



duction see Brown et al. (2007). 

Two corrections have been applied to our sam- 
ple of extragalactic sources. Firstly, structure 
functions are sensitive to large-scale gradients in 
electron density across the field of view and in- 
clude a geometrical component due to the change 
in direction of the regular magnetic field, (see e.g. 



Brown et al.ll2007h . As a first order correction, 



we approximate this with a 2D linear gradient in 
RM, and subtract this from the region over which 
a structure function is computed. Furthermore, 
lines of sight through discrete structures like H II 
regions and supernova remnants can have decep- 
tively large RMs due to an increased electron den- 
sity and possibl y magnetic field wi thin these local- 



ized structures ( Mitra et al. I l2003h . Therefore, we 



have used the total intensity 1.4 GHz radio data 
from the ATCA combine d with Parkes sin gle-dish 
data as well as Ha maps (lFinkbeinerll2003h to dis- 
card 27 extragalactic sources with a sight line pass- 
ing through a visible supernova remnant or H II 
region. These sources are listed in Table [T] We 
recognize that omitting RMs through discrete high 
density regions may introduce a bias in the struc- 
ture function, viz. decrease the structure function 
amplitude on large scales. However, the results are 
very similar to results without discarding extreme 



RMs (cf. I Haverkorn et al. 2006a ). indicating that 
the bias, if present, is low. 

The SGPS data probe the inner Galaxy, which 
includes a number of spiral arms. Consequently, 
these data are well-suited to study differences 
in the structure in the ISM in spiral arms and 
in intcrarm regions. We constructed second- 
order structure functions of RM, Drm(59) = 
<[RM(6») - RM(6 + S0)] 2 ) g , for different lines of 
sight in 'spiral arms', i.e. lines of sight primarily 
through spiral arms, and 'interarm regions', sight 
lines mostly through interarm re gions, estimated 
from the spiral arm positions in Cordes fc Laziol 
(|2002l ). The lines of sight used to separate spiral 
arms and interarm regions are shown in Fig. [1] 
The error bars denote errors propagated from un- 
certainties in the RM values. For lines of sight 
at high longitudes close to the Galactic center it 
is not possible to distinguish between 'spiral arm 
lines of sight' and 'interarm lines of sight', because 
spiral arms start running perpendicular to the line 
of sight. Therefore, we do not use data at I > 326° 
in this analysis. 

Figure [2] shows Drm for spiral arms and in- 
terarm regions. The fi gure is similar to Fig. 1 in 
Haverkorn et al. ( 20061 ) except here we have used 



3 http: / /www. atnf.csiro.au/research/cont/sgps/queryForm. html 



the finalized RM list in iBrown et alj (|2Q07I ) , and 
discarded the 27 sources in Table Q] The spiral 
arm structure functions are flat, while in interarm 
regions the structure functions rise, and in two out 
of three cases show a turnover from a power law 
at small scales to flat at the larger scales. The lo- 
cation of the turnover is interpreted as the largest 
angular scale of structure in the interarm regions. 
Using the argument that the largest angular scales 
in RM are probably coming from nearby, assum- 
ing a distance of 2 kpc yields an outer scale at a 
spatial scale of about 100 pc. For the spiral arms 
we can only say that the outer scale of structure, 
i.e. the smal lest scale we probe, is smaller than 
about 10 pc (jHaverkorn et al.ll2006a[) . 

The amount of sources used to estimate the 
structure functions are 50, 20 and 18 for the re- 
spective interarm regions, and 8 sources for each 
of the spiral arms. Given the low source density 
in the spiral arms, we assess the reliability of the 
results here. The bin size is somewhat restricted 
due to the paucity of sources, but bin sizes between 
0.5° and 1° are reasonable and yield comparable 
results. 
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Table 1: Extragalactic point sources discarded from the analysis because they located are behind H n regions 
or supernova remnants, detected in total intensity or Ha. The columns give the longitude I and latitude b 
of the source in degrees, its RM in rad m~ 2 and the reason for flagging it: Stokes I or Ha. 



Figure [3] shows the structure functions of RMs 
towards the Carina and Crux arms without any 
binning of sources. The solid lines give linear fits 
of the data, confirming our statement that the 
structure function in the spiral arms is flat. In 
the Carina arm, the seven uppermost points (all 
at log(£>AA/) > 5.3) are the combination of one 
source with extreme RM and all other sources in 
the region. The presence of this extreme source 
makes clear why the two data points on the largest 
scales are lower than the other points: on these 
scales the extreme RM source does not contribute. 
Furthermore, this explains why the amplitude of 
the structure function is higher in the Carina arm 
than in the Crux arm: omitting this source yields 
comparable amplitudes for both arms. 

However, we are hesitant to discard sources on 
the basis of extreme RM alone, as there is no rea- 
son per se why these sources would not be part 
of the spectrum. Therefore, we only omit sources 
visibly behind a discrete structure, as discussed 
above, and leave this extreme RM source in the 
dataset, while commenting on changes when this 
source is omitted. 

We can estimate the outer scale from modeling 
of the amplitude and slope of the structure func- 
tions, or from the amount of depolarization of the 
point sources by the Galaxy. These two methods 
will be discussed in Sections [3] and |H respectively. 



3. Outer scale from rotation measure 
structure functions 

The structure function slopes in Fig. [2] are 
0.32±0.05, 1.09±0.18 and 0.71 ±0.17 for interarm 
regions 1, 2, and 3, respectively, much shallower 
than those expected from a Kolmogorov structure 
function, which has a slope m — 5/3. However, we 
argue in this section that on smaller scales the RM 
structure function has to turn over to a steeper 
slope to be consistent with observations of elec- 
tron density fluctuatio ns on smaller scales (e.g., 
Armstrong et al.l ll995); additional arguments for 
a steeper RM power spectrum on smaller scales 
are given in Section [5l 

Minter & Spangler (1996, hereafter MS96) de- 
veloped a formalism to describe the structure func- 
tion of RM assuming power spectra in magnetic 
field fluctuations and in electron density fluctu- 
ations which are zero-mean, isotropic, Gaussian, 
and independent. Assuming Kolmogorov turbu- 
lence the RM structure function, Drmi can be 
described as: 
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Fig. 1. — Bird's-eye view of the Galaxy, with the 
longitude ranges that define the 'spiral arms' and 
'interarm regions'. The distance L which contains 
90% of the total electron density along the line of 
sight is given by the length of the dashed line. The 
distances 0.5L at which the regular magnetic field 
strength is calculated is given by the black crosses. 
The dotted circle denotes the distance L* used for 
calculating angular outer scales. 
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(1) 



where n e0 is the mean electron density, B || is 
the mean magnetic field strength along the line 
of sight, Zq the outer scale of the Kolmogorov tur- 
bulence, and L the length of the line of sight. The 
coefficients C% and are defined in the descrip- 
tion of the magnetic field and density fluctuations 
as power laws with the same outer scale Iq and 
spectral index a such that 

(^(r )^(r + r))= J d 3 q-^^ (2) 

where wave number q — 2ir/r and go = 2tt/Iq. A 
similar expression applies for (5n(r )Sn(r +r)). 
The spectral index of the power spectrum is a = 
11/3 for Kolmogorov turbulence, and is related to 
the slope of the structure functio n Drm oc r m 
asr7i = a- 2for2<a<4 (iRickettl [l977l) . 



Fig. 2. — Structure functions of RM for Galactic 
interarm regions (top) and spiral arms (bottom). 
The solid lines are linear fits to the rising parts of 
the structure functions. 



is only true on scales smaller than what we ob- 
serve. However, arguing that the slope turns over 
to Kolmogorov on smaller scales, we use a Kol- 
mogorov dependence according to equation (TTJ) on 
scales I < Iq , which flattens to the shallow or flat 
observed spectra for I > Iq . 

Most of the input parameters in equation (JTJ 
are known, so that by joining this equation 
with RM structure functions in Figure [2] the 
outer scale of Kolmogorov turbulence Iq can 
be computed. The amplitude of the electron 
density fluctu ations C~j is ta ken to be C"i = 



io- 



.-20/3 



([Armstrong et alJ 119951). and the 



Equation [T] is only valid if 89L/Iq < 1, which 



magnetic field fluctuations C% can be derived 
from the observed value from MS96 as C% = 
5.2B2 an (^)-2/3 /iG 2 m -2/3 ; where B ran is the 

strength of the random component of B in /iG 
and Iq is in parsecs. Values for the mean elec- 
tron density n e Q, mean magnetic field B , random 
magnetic field B ran , and distance to the emission 
will be derived in the following subsections. 

3.1. Magnetic field 

A number of estimates of the total magnetic 
field strength in the Galaxy, based on a number 
of observations, indicate that the total magnetic 
field strength is around 6 nG at the Solar radius. 
Heilesl (Il995h gives an extended discussion about 
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Fig. 3. — Structure functions of RM for the Carina 
spiral arm (left and the Crux arm (right). Each 
data point denotes a single source pair. The solid 
lines denote linear fits of all data. The dashed 
line in the Carina arm plot indicates a linear fit 
after discarding of the upper row of points at 
log(D RM ) > 5.3, which are all caused by one 
source of extreme RM. 



the different ways to determine Galactic magnetic 
field strengths, i.e. using the synchrotron emissiv- 
ity under the assumption of minimum energy or 
minimum pressure, or with the cosmic ray den- 
sity measured in the solar neighborhood, and us- 
ing Faraday rotation from pulsars. His estimates 
for the total magnetic field in the solar neighbor- 
hood range from 4 to 7.4 fiG for different assump- 
tions, while he estimates B to t ~ 7.6 — 11.2 /zG 
at Gala ctocentric rad i us Rg „i — 4 kpc. More 
recently. IStrong et al. ( 2000T ) modeled cosmic ray 
evolution and propagation in the Milky Way, us- 
ing constraints from synchrotron and 7-ray emis- 
sion. Their results indicate a total magnetic field 
strength of 6.1 fiG at the solar circle, increasing 
exponentially towards the inner Galaxy with a 
scale length of 10 kpc. These measurements are 
in good agreement with ea ch other, with re sults 
from synchrot ron radiati on jgeck et al. 1996) and 
from pulsars (jHan et al.ll2004h . 

We can estimate the relative strengths of the 
magnetic field in the s piral arms and interarm 
reasons from Fig. 7b in Beuermann et all (|l985h ■ 
which shows that the mean volume emissivity in 
the spiral arms is twice as high as the mean volume 
emissivity in the interarm regions, i.e. e arms = 
^interarms = 1.5(e), where the average volume 



arms must be (3/2) 2 / 7 that in the interarm re- 
gional where the exponent 2/7 occurs because in 
the minimum energy approximation (B) oc e 2 / 7 . 
We adopt the dependence of the magnetic field 
streng th on Galactocentric radius in IStrong et al 



(20(KJ), and correct for the relative strengths in the 
arms and interarm regions as given above, lead- 
ing to the total field strengths given in Table [2] 
The regul ar magneti c field component B reg is esti- 



mated bv lHan et al. ( 20061 ) from pulsar dispersion 



and rotation measures as B reg = 2.1 ± 0.3 /xG, in- 
creasing exponentially inwards with a scale length 
of 8.5 kpc. We adopt these values for the reg- 
ular magnetic field strength, so that B ran = 

J Bf ot — B^ eg w 0.94£? tot , independent of Galac- 
tocentric radius. The component of the regular 
magnetic field parallel to the line of sight _B || in 
equation (pj is evaluated through B || = -B cos/3, 
where (3 is the angle between the line of sight and 
the re gular magne tic field B . A pitch angle of 
— 12° (| Vailed 120041 ) is assumed, but a pitch angle 
of 0° gives no significant changes. 

3.2. Electron density 

The average electron density used for each line 
of sight is d erived from the NE200 1 electron den- 
sity model ( Cordes fc Laziol 120021 ) . The electron 
density was evaluated from the model for a par- 
ticular line of sight centrally through each arm 
and each interarm region. The adopted average 
electron density is the average density over the 
adopted line of sight. Most of the structure in 
NE2001 is on such large scales that a change in 
direction of the line of sight of a degree or so does 
not influence the results significantly in any of the 
lines of sight. The adopted values for the mean 
electron density along each line of sight can be 
found in Tabled 

3.3. Distances 

Two different distances for modeling are needed: 
L, the total path length of the Faraday rotating 
material, and L*, the distance used to convert the 



emissivity is (e) = 0.5(e arms + e interarms ). There- 
fore, the total magnetic field strength in the spiral 



9 We follow the argument put forth in[Heilcs ( 1995), although 
he uses Earms ~ 2(e) = 2(e arms ~\~ Winter arms) ■ His approx- 
imation tjnter arms ~ is reasonable i n the outer Galaxy. 
However, the iBeuermann et all J1985T) synchrotron maps 
indicate that e ar ms = ^interarms is more appropriate in 
the inner Galaxy. 
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Fig. 4. — Estimate of the structure function of 
RM, Drm, for a constant magnetic field (i.e. 
Bran = 0), with the extrapolation from electron 
density fluctuations measured on smaller scales for 
the input parameters of the Carina and Crux arms 
(dotted lines). The asterisks are data from the 
Carina and Crux arms, with solid line fits. The 
dashed lines show the suggested connection to the 
data and the arrows denote the outer scale iff . 



outer scale of fluctuations from angular to spatial 
scales, both of which are shown in Fig. [TJ The 
distance L is chosen as the distance to the point 
for which 90% of the electron density is contained 
along the path, and is given in Table [U Subse- 
quently, the distance from the sun at which B to t (r) 
is computed for each sight line is 0.5L. 

If the statistical properties of the medium do 
not change along a given line of sight, the distance 
at which a certain angular scale is expected to cor- 
respond to the largest spatial scale is a small dis- 
tance. Consequently, L* , which corresponds to the 
largest spatial scale iff , should be different from L, 
and is estimated to be 2 kpc. 

These distances are very rough estimates. How- 
ever,we show in Section 13.51 that the sensitivity 
of our conclusions to the anticipated error in dis- 
tances is low. 

3.4. Results from RM structure functions 

We extrapolate the observed slopes to smaller 
scales, with a steepening to a Kolmogorov spec- 
trum at scale iff , which is the outer scale of the 



Kolmogorov part of the spectrum. The constraint 
of equal amplitudes between the Kolmogorov 
structure function on the small scales and the 
shallower structure function on the larger scales 
at turnover scale iff yields 



Drm(60) 



505/3 

(lff/L*) 5 / 3 ~ r 



for 
for 



< lff/L* 
> lff/L* 



(3) 

where m is the spectral index of the shallower 
structure function. In reality the structure func- 
tion will not make a sharp break as described here 
and instead will sho w a gradual turnover o ver a 



range of scales (e.g., lHaverkorn et al.l 12004 ). but 



we use this parameterization as a good first ap- 
proximation. Combining equations |T]) and J3|) for 
89 > lff/L* yields 



D RM (S9) 



[ Cin 2 e0 Bt ot 
+C 3 B to t cos j3' 



K\2/3 



xi 8/3 ( *0_ 

L* 



5/3— m 



SO" 



(4) 



where d = 488.37, C 2 = 0.04 and C 3 = 0.08 
are constants derived from the known variables 
in equation |T]), and m is the spectral index of 
the structure function. We fit equation (j4|) to the 
observed structure functions to obtain estimates 
for the outer scale of the Kolmogorov part of the 
power spectrum ^(SF), with results given in Ta- 
bled 

To compare our results with the Kolmogorov 
spectra found in electron density, we evaluate the 
structure function of RM caused by electron den- 
sity fluctuations only (B ran — 0), shown as the 
dotted lines in Fig. [5] for the input parameters of 
the Carina and the Crux arm. This is a lower limit 
for the structure function of RM if the power spec- 
trum of electron density on small scales is extrap- 
olated to the large (parsec) scales discussed here. 
The spiral arm data, plotted as asterisks in Fig. [U 
fall below the extrapolation of D rm to large scales, 
indicating that Drm has to turn over at smaller 
scales where the dotted and dashed lines in Fig. 2] 
meet. The same argument holds for the interarm 
data, which are not shown in the Figure for clar- 
ity. Although the input parameters are uncertain, 
this result is remarkably stable against variations 
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of the input parameters (Section 13. 5p . For fluctu- 
ations in electron density only, this turnover scale 
is around 3 pc (1 pc) for the Carina (Crux) arm. If 
magnetic field fluctuations are present as well, the 
turn over scale Iq is a little smaller, as noted in 
Table [2] Therefore, the Kolmogorov spectrum in 
electron density on small scales, if extrapolated to- 
wards larger scales, does not extend all the way up 
to scales of ~ 100 pc as previously assumed, but 
displays a break to a shallower slope (interarms) 
or a constant value (spiral arms). If it did, equa- 
tion ([4} demonstrates that at a scale r — 100 pc, 
Drm would be 2.3 orders of magnitude higher 
than the observed values. 

3.5. Sensitivity of results to input param- 
eters 

Due to the non-straightforward dependence of 
Drm on Iq in equations and (01, we tested 
numerically how sensitive the results are to varia- 
tions in the input parameters L, L* and n e . If 
the path length L is increased or decreased by 
30%, this will decrease or increase the resulting 
outer scale Iq by a factor of 50% . The same ef- 
fect is seen for an increase or decrease in L* or 
n e by a factor two. Although it makes sense to 
assume that the largest angular scales S8 cor- 
respond to the largest spatial scales Iq at some 
nearby position along the line of sight, as we have 
done, even if we assume that the largest spatial 
scales are at the midway distance along the line 
of sight, Iq- = 0.5 Ltan(<56>o), the obtained outer 
scale Iq » 1 — 5 pc for both arms and interarm 
regions. Also, if the amplitude of electron density 
C% were a factor 10 differ e nt in our data from the 
data in lArmstrong et al. (1995). the outer scale 
would change by less than 50%. As the results are 
fairly robust against reasonable changes in the in- 
put parameters, we feel confident in asserting that 
the outer scale of Kolmogorov turbulence Iq must 
be on the order of a few parsecs. We note, how- 
ever, that due to the assumptions made and un- 
certainties in input parameters, a relatively large 
uncertainty in the determined outer scale has to 
be taken into account. 



4. Outer scale from depolarization of 
point sources 

An independent estimate of the outer scales 
of magneto-ionic structure can be made from 
depolarization of extragalactic point sources, if 
caused by structure in the foreground ISM on 
angular scales smaller than the size of the (un- 
resolved) source. Intrinsic variations in polariza- 
tion angle causing partial depolarization are ex- 
pected to arise within any polarized extragalactic 
source. Indeed, no source in our sample exhibits 
the theoretical maximum degree of polarization 
of around 70% (|Pacholczvklll970l ) but instead ob- 
served degrees of polarization lie typically below 
10%. Depolarization by foreground components 
can be caused by beam depolarization due to mag- 
netic field and/or electron density fluctuations 
on scales smaller th an the source size (see e.g. 
Gaensler et al.ll2005) or by bandwid th depolariza- 



tion (see e.g. IStil fc Taylorl |2007[) . Bandwidth 



depolarization is given by p = pO sin(Afl)/Afl 



wher e A0 = 2RMc 2 Av/v 3 ([Gardner & Whiteoak 
1966). However, significant bandwidth depolar- 
ization across our frequency band Au (8 MHz) 
can only be achieved by RM > 5700 rad m~ 2 , 
much higher than observed RMs. Therefore, any 
foreground depolarization in our data is caused by 
beam depolarization across the face of the source 
rather than bandwidth depolarization. 

Figure [5] shows RMs in the upper panel and 
degree of polarization, p, in the bottom panel. A 
clear anti-correlation between |RM| and p is visible 
especially at the lower longitudes. As the scale of 
the structure in RM and p is several degrees, this 
cannot be intrinsic to the sources but instead must 
be caused by the Galactic ISM. T his agrees with 
the results of lGaensler et aL ( 2005 ) who found an 
anticorrelation between degree of polarization and 
Ha, which is correlated with RM|. Due to the 
power law nature of RM as a function of scale, a 
high RM also indicates large fluctuations in RM. 
This is shown in Figure [6] which shows the stan- 
dard deviation in RM as a function of fractional 
polarization, for all data together and for the spi- 
ral arms and interarm regions separately. As ex- 
pected, the data in the spiral arms show a higher 
standard deviation in RM and a lower fractional 
polarization than in the interarm regions. 

The Galactic component of depolarization can 
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Fig. 5. — The SGPS field of view, where circles denote the position of extragalactic sources. Top: Filled 
(open) symbols represent positive (negative) RM scaled linearly as —1000 rad m~ 2 < RM < 1000 rad m~ 2 . 
Bottom: symbols represent the degree of polarization, p, scaled between 0.4% and 13.7%. 



be estimated from the power law spectrum of RM 
fluctuations. The angular size of the outer scale 
of structure 8q is much larger than the angular 
source size 6 src . In this approximation 6q >> 
9 S rc, the depolarization by a power spectrum of 
RM fluctuations is give n by th e degr ee of polar- 
ization p, adapted from iTribbld (Il99lh as: 



<(^) 2 ) » 1 - 4cx 2 A 4 2"V 2 ( ^Y" F(l + 
Po V Iq J 2 

(5) 

where po is the degree of polarization of the extra- 
galactic source when its radiation arrives at the 
Milky Way, and r src is the size of the source. The 
RM standard deviation a, structure function slope 
m and outer scale Iq are related via the structure 
function D RM {r) = 2a 2 (r/l$) m for r < lg . 

A decrease in degree of polarization p with in- 
creasing standard deviation in RM er, as predicted 
by equation [5l is visible in Figure [6] for all data 
together and for the spiral arms and interarm re- 
gions separately. As expected, the data in the spi- 
ral arms show a higher standard deviation in RM 
and a lower fractional polarization than in the in- 
terarm regions. With an estimate of po and r src , 
we can fit a as a function of p to the data points, 
and obtain a best-fit outer scale Iq . 

The amount of intrinsic depolarization result- 
ing in polarization degree po can be estimated 
from the extragalactic sources observed around 
the LMC (Gaensler et al. 2005) to be 10.4%, which 



we assume is the average degree of polarization 
of unresolved point sources at 1.4 GHz for which 
all depolarization is intrinsic. With these assump- 
tions, the depolarization below 10.4% is then due 
to the variations in Galactic RM across the face of 
the source with a median size of 6 arcsec for this 
flux density range (Gaensler et al. 2005). Note 
that the percentage of intrinsic depolarization is 
higher than the actual average degree of polariza- 
tion due to a selection of strong, highly-polarized 
sources over weak, weakly polarized ones, and due 
to selection of sources with linear </>(A 2 ) behavior. 
However, since we are interested in the relative 
depolarization only, this selection effect does not 
influence our conclusions. 

In the spiral arms it is straightforward to use 
equation (JSJ) to determine the outer scale Iq 
needed to obtain the observed depolarization. As- 
suming Kolmogorov turbulence (m = 5/3), we 
determine the value of the RM standard deviation 
a from the structure function saturation level. For 
the interarm regions we observe a spectrum which 
is considerably shallower than a Kolmogorov spec- 
trum. Assuming that this spectrum turns over 
to a steeper Kolmogorov spectrum towards small 
scales, the Kolmogorov slope on small scales will 
dominate the depolarization of the point sources. 
Therefore a can be assumed to be the value of the 
structure function at the scale Iq , which is the 
outer scale of the Kolmogorov turbulence, i.e. the 
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scale at which the Kolmogorov slope turns over 
into a shallower slope. 

The data are best represented by equation ([5|) 
for iff = 0.2°, 0.25°, and 0.1° for all data, spiral 
arms and interarm regions, respectively (shown in 
Figure H as a solid line, with outer scales twice 
higher and lower indicated by dashed lines), so 
that the outer scales of Kolmogorov turbulence 
are approximately 8.7 pc in the spiral arms and 
3.5 pc in the interarm regions. The data in the spi- 
ral arms are fairly sparse - confirming additional 
depolarization in the spiral arms - so that the stan- 
dard deviations and the fit to the model are un- 
certain. However, in the interarm regions the de- 
polarization model is a good fit to the data, and 
the probability that the standard deviations are 
constant with fractional polarization is < 0.1%. 

The estimates of iff from depolarization are 
somewhat larger than iff from the structure func- 
tion analysis. However, if we consider that the 
errors in the distances are large, in conjunction 
with assumptions such as the lack of correlation 
between magnetic field and electron density, the 
difference in estimates of iff from the two meth- 
ods is not necessarily significant. Certainly both 
methods indicate that the outer scale of the Kol- 
mogorov spectrum is likely a few parsecs, much 
smaller than the previously assumed value of ~ 
100 pc. 

5. Expected steep spectrum on small 
scales 

As argued in the previous section, the RM 
structure function has to turn over to steeper 
slopes on smaller scales to be consistent with the 
electron density fluctuation data on smaller scales. 
There are two additional reasons why the RM 
structure functions cannot continue to have the 
same shallow or flat slope on smaller scales. 

Firstly, a steep Kolmogorov-like magnetic field 
power spectrum is indicated by cosmic ray data. 
As cosmic rays are most effectively scattered by 
magnetic field fluctuations on the same scale as 
their ion gyro radius, cosmic ray losses as a func- 
tion of energy are closely related to the magnetic 
field power spectrum. The cosmic ray distribu- 
tion as calculated from a leaky box model can ex- 
plain cosmic ray observational data if the Galac- 
tic magnetic field has a power spectrum with a 
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fractional polarization p 

Fig. 6. — Diamonds show standard deviation of 
RM, cr, versus degree of polarization, p, for all 
data (top), spiral arms only (middle) and inter- 
arm regions only (bottom) . Sources are binned in 
degree of polarization with bin width Ap = 0.02. 
The solid lines are model predictions from equa- 
tion |(5]), for iff = 7.7 pc, 8.7 pc and 3.5 pc for 
the top, middle, and bottom panels, respectively. 
The dashed lines show the models for iff a factor 
2 larger and smaller. 



Kolmogorov spectral index (jJokipii 1988). Fur- 
thermore, the cosmic ray power spectrum is re- 
markably smooth^! on scales of 10 9 to 10 18 eV, 
corresponding to gyro radii of 10 20 to 10 12 cm. 
Therefore, the magnetic field power spectrum is 
also expected to be smooth on these scales. 

However, recent numerical simulations show 
that the magnetic field power spectrum does not 
necessarily follow the electron density spectrum. 
In fact, magnetohydrodynamic simulations in the 
limit of a weak homogeneous magnetic field show 
that the magnetic field fluctuations are all concen- 



10 The steepening of the spectrum around 3 X 10 15 eV (the 
"knee") and flattening around 3 X 10 18 eV (the "ankle") 
are very slight and not relevant for this argument. 
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trated on scales much smaller than thos e under 
discussion here ( Schekochihin et al. 2002^ 1 . This 
would indicate a flat structure function of mag- 
netic field on larger scales. At first sight, this the- 
ory agrees with the observations shown in Fig. [2J 
where the weak-mean-ficld approximation could 
be applicable to the spiral arms, and structure 
functions would be expected to be flat. How- 
ever, this would indicate that magnetic field fluc- 
tuations of a few microgauss would be present on 
scales as small as a fraction of a parsec. In this 
case the degree of polarization, p, as a function of 
the intrinsic degree of polarization, po, is given by 

P = Po e-^ 4 *e-™ (6) 

(|Burnl Il966t ISokoloff et~ai] Il998t ) in the approxi- 
mation that the scale of the fluctuations is much 
smaller than the telescope beam. The standard 
deviation of RM, a, is derived from the flat struc- 
ture functions in the spiral arms. Therefore, mag- 
netic field fluctuations on these scales would com- 
pletely depolarize the synchrotron radiation in 
the Galaxy at 1.4 GHz. Since we observe polar- 
ized radiation at this frequen cy coming from all 
over the Galactic pla ne (e.g. iTavlor et al.l [2003; 
Haverkorn et al.ll2006l ) , magnetic field fluctuations 
cannot remain at this magnitude towards smaller 
scales. Instead, we showed in Section [4] that the 
observed amount of depolarization is consistent 
with a RM structure function with a Kolmogorov 
slope and an outer scale of a few parsecs or smaller. 

6. Discussion 

The conclusion that the outer scale of turbu- 
lence in the spiral arms is observed to be on the 
order of a few parsecs is not expected or straight- 
forward. Based on evidence discussed in the intro- 
duction, outer scales of turbulence are expected to 
be on the order of 100 pc, which is one or two or- 
ders of magnitude higher than those observed here. 
Earlier work noted the small outer scale of fluctua- 
tions in the spiral arms, but attributed that outer 
scale likely to H II regio ns along the line of sight 
( Haverkorn et al. 2006a[ ). H II regions are ubiqui- 
tous enough in the spiral arms to dominate fluctu- 
ations in RM, and this solution would reconcile a 
larger outer scale of turbulence with a smaller ob- 
served outer scale of RM fluctuations. However, if 
this were the case, the amplitude of the RM struc- 
ture functions would lie above the amplitude of the 



turbulence structure function. Our estimate of the 
amplitude of the RM structure functions indicates 
that the observed structure functions lie far below 
the lower limit for the RM structure function if 
extrapolated from electron density fluctuations on 
small scales. If the fluctuations in RM associated 
with Kolmogorov turbulence were to continue up 
to scales of a hundred parsecs, the observed am- 
plitude of the structure function of RM would be 
much higher than that observed. 

In the interarm regions, a plausible option is 
multiple scales of energy input: for supernova- 
driven turbulence, the outer scale is believed to 
be about 100 pc (as observed). However, if en- 
ergy sources such as stellar winds or outflows, 
interstellar shocks or H II regions input a sig- 
nificant amount of energy into the interstellar 
turbulence on smal l er sca les (typically parsecs, 
Mac Low fc Klessenl (|2004l )). this may flatten the 
structure function on scales of ~ 1 pc to scales of 
~ 100 pc, consistent with our observations. 

This does not contradict earlier studies that 
reported outer scal e s of t he order of 100 pc. 
Simonetti fc Cordesl ( 1986h found larger outer 
scales, but their data included large parts of the 
sky, most or all of wh i ch we re located at higher 
latitudes. IClegg et al. (|l992h present some data 
in the Galactic plane, but the outer scale of fluc- 
tuations in those data is not well determined due 
to the paucity of data and the included geomet- 
rical component of the magnetic field. Leahy's 
(1987) structure functions of RMs of sources in 
the Galactic plane are consistent with an outer 
scale of a parsec, as is the analysis of DM vari- 
ations of close pulsar pai r s in globular clusters 
(jSmirnova k. Shishov|[20o3 lRansomll2007l ). 

So the picture arises of a smaller energy input 
scale of turbulence or fluctuations in the magneto- 
ionized ISM in the Galactic plane, while larger- 
scale structure exists in the Galactic thick disk or 
halo. Indications of increasing correlation lengths 
of the magnetic field with height above the galaxy 
case hav e been found in a number of external 
galaxies (DumkeetaL I ll99ffl . 



7. Summary and conclusions 

Faraday rotation measurements of polarized ex- 
tragalactic sources behind the inner Galactic plane 
have been used to study the characteristics of 
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long 


;itude 


Btot 


n e o 


L 


L* 




iff (depol) 




1 


:°] 


\MG] 


[cm" 3 ] 


[kpc] 


[kpc] 


M 


[pc] 


Interarm 1 


255 


- 280 


4.5 


0.03 


11.5 


2 


2.3 


3.5 


Interarm 2 


290 


- 305 


4.6 


0.045 


16 


2 


0.8 


8.7 


Interarm 3 


315 


- 326 


5.2 


0.075 


18 


2 


0.3 


3.5 


Carina arm 


280 


- 290 


8.3 


0.06 


14.5 


2 


2.4* 


8.7 


Crux arm 


305 


- 315 


9.0 


0.07 


17 


2 


1.0 


3.5 



* This value is 0.65 pc if the extreme RM source discussed in Section [2] is omitted. 



Table 2: ISM parameters for three interarm regions and the Carina and Crux spiral arms. The parameter 
n e o is the electron density averaged over the line of sight; L and L* are lines of sight as described in the 
text; and B to t is the total magnetic field strength. The outer scale of the Kolmogorov power spectrum iff is 
given as obtained from structure functions (SF) and from depolarization of point sources (depol). 



the magnetized, ionized interstellar medium in the 
plane, in particular in the spiral arms and in in- 
terarm regions. Rotation measure structure func- 
tions show a shallow slope in the interarm regions 
and saturation on a scale of ~ 100 pc, i.e. there 
are no fluctuations on scales larger than the satu- 
ration scale. Flat structure functions in the spiral 
arms indicate that the outer scale of RM fluctua- 
tions in the spiral arms is smaller than ~ 10 pc, 
the smallest scale observed. These shallow and 
flat structure functions must turn over to steeper 
slopes towards smaller scales for three reasons: (1) 
to match up with the electron density power spec- 
trum on subparsec scales, assuming the large and 
small scale datasets are part of the same power 
spectrum; (2) a shallow RM structure function on 
smaller scales would give more depolarization than 
observed; and (3) cosmic ray distribution data and 
the smooth cosmic ray power spectrum indicate 
a smooth magnetic field power spectrum with a 
slope similar to the Kolmogorov slope. The scale 
of the break in the structure function is the outer 
scale of the Kolmogorov power spectrum iff , and 
is estimated using two independent methods: the 
analysis of RM structure functions, and by model- 
ing the depolarization of the extragalactic sources. 
Given the large uncertainties in input parameters, 
both methods agree reasonably well and imply an 
outer scale of the Kolmogorov slope of a few par- 
sees. This estimate is almost two orders of mag- 
nitude smaller than the generally assumed outer 
scale of ISM turbulence of ~ 100 pc. However, 
extrapolating the observed electron density fluc- 
tuations on small scales to parsec scales shows 
that the amplitude of the structure function would 



be orders of magnitude higher if the Kolmogorov 
spectrum did in fact extend out to 100 pc. Instead, 
the outer scale of Kolmogorov turbulence iff that 
we obtained from our observations indicates that 
energy in ISM turbulence is injected on scales of a 
parsec rather than 100 pc. This is the main energy 
injection scale in the spiral arms, which show flat 
structure functions on scales larger than that. In 
the interarm regions, the structure functions keep 
rising although not as steep as Kolmogorov turbu- 
lence, indicating an additional source of structure. 
We propose that in the spiral arms stellar energy 
sources such as stellar winds and protostellar out- 
flows are the predominant sources of turbulence, 
whereas in the interarm regions there is evidence 
of energy injection on larger scales, most likely 
caused by supernova remnant and superbubble ex- 
pansion. 
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